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We examine the dynamical behavior of matter coupled to gravity in the context of a linear 
Klein-Gordon equation coupled to a Friedman-Robertson- Walker metric. The resulting ordinary 
differential equations can be decoupled, the effect of gravity being traced in rendering the equation 
for the scalar field nonlinear. We obtain regular (in the massless case) and asymptotic (in the 
massive case) solutions for the resulting matter field and discuss their ensuing finite time blowup in 
the light of earlier findings. Finally, some potentially interesting connections of these blowups with 
features of focusing in the theory of nonlinear partial differential equations are outlined, suggesting 
the potential relevance of a nonlinear theory of quantum cosmology. 

In the past two decades there has been an extensive interest in 3-dimensional gravity, particular after the demon- 
stration of the fact that its quantum version is solvable [1] and that it contains black hole solutions [2]. 

The idea of examining cosmological models in this context was, however, of interest before [3-5], as well as after 
[6,7] these findings. In most of these studies [3-5,7], the Friedman- Robertson- Walker (FRW) metric was used and the 
resulting equations were ordinary differential equation governing the time-evolution of the scale factor of the relevant 
OO ' metric and the evolution of the matter and/or radiation field coupled to it. 

On a slightly different track, one can list the works of [6,8,9] (see also references therein and the review of [10]), where 
the metric scale factors were allowed to be temporally as well as spatially variable, and the resulting partial-differential 
^ , equations (PDEs) were studied to obtain collapse type solutions. 

In this brief report, we will restrict ourselves to the former type of considerations (but in a 3+1 dimensional setting), 
namely fixing the spatial dependence of the metric tensor to be of the FRW type (in co-moving coordinates and with 
,— ' ' a "cosmological" time choice) 
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ds 2 = ~dt 2 + a 2 (t) ( 3_dr 2 +r 2 d0 2 + r 2 sin 2 Odd) 2 ] (1) 
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O and coupling gravity to the simplest possible model for matter, namely a linear Klein-Gordon (KG) type equation for 
a massive scalar. In the metric of Eq. (1), a is the scale factor while k describes the curvature of the spatial slice and 
5h i can be normalized to the values — 1, 0, 1 in the hyperbolic, flat and elliptic case respectively. Previously such models 
^P" were examined in a more complicated setting (most often in 2+1 dimensions) where either equations of state [7] or 
, special scalar field potentials [4] were used to obtain closed form solutions. 

Here, we will examine the simplest possible (physically relevant) scalar field potential V((f>) = m 2 (f> 2 /2, which gives 
rise to a linear (in </>) equation for the scalar coupled to gravity. In particular the Einstein-Klein-Gordon field equations 



X 

in this case read: 



k + a 2 

-3^^+P0 = O 2 
a 2 + 2ad 



p = (3) 

4> + m 2 (f)+—(f) = (4) 
a 

where the energy density and the pressure associated with the scalar field are given respectively by (see e.g., [7]) 

P^\{4> 2 +m 2 4> 2 ) = \4> 2 + V{cp), (5) 
p<> = \{4> 2 -m 2 4> 2 )^\tf -V{4>). (6) 
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Note that the first is the quadratic constraint Gg = Tq ; the second is the only independent spatial equation G\ = T± 
while the last is the dynamical equation for the scalar field, i.e., the Klein-Gordon equation. G and T denote the 
Einstein curvature tensor and the energy-momentum tensor respectively. 

In the above equations the dot denotes temporal derivative. Naturally, Eqs. (2)-(4) are not linearly independent 
as the linear combination of the derivative of the first and of the second can be used to obtain the third. 

We first examine the spatially flat case of k = 0, and use it as a guide. In this setting one immediately observes 
that from Eqs. (2) and (4), two separate expressions for a/ a can be obtained, hence equating them, a second-order 
ordinary differential equation (ODE) emerges for the scalar field <fr in the form: 

^(mV+^(m^^) 2 . (7) 

Notice that this is the only case (among the ones that we will examine) in which the resulting ODE is of 2nd order. 
In the remaining cases, the ODE is of 3rd order. Moreover, it is worth commenting on the nature of this ODE: 
in particular, the resulting equation is nonlinear. Hence, even though the inclusion of matter is realized through a 
linear dynamical equation, in the setting of even the simplest cosmological models, the coupling of matter to gravity 
induces the emergence of a nonlinear equation; it is as if the "trace" of gravity, when the equation for the scalar field 
is decoupled from it, remains in the nonlinearity of the resulting ODE. 

As the simplest possible case among the ones with k = 0, we examine the massless case, i.e., m = 0. In the latter 
setting, we obtain the solutions (up to a constant shift) of the form: 

= ±y|log(t*-t) (8) 

which in turn results in a power law dependence of the scale factor of the form a <~ (t* — i) 1 / 3 . Hence, a blowup 
(focusing) type effect occurs at t = t* in a logarithmic fashion for the scalar and the corresponding scale factor shrinks 
to with a power law dependence. Such dependencies are reminiscent of the critical blowup in prototypical nonlinear 
PDEs with focusing solutions (for a spatio-temporally dependent field tp) such as the nonlinear Schrodingcr equation 
[11-13] 

ii/> t = -AV> - M 2ct ^. (9) 

A stands for the Laplacian and a the power of the nonlinearity. For da < 2, no focusing solutions occur (d is 
the dimensionality of the Laplacian); when da = 2, logarithmically slow focusing phenomena take place (in the 
corresponding "proper time" see [11-13]), while when da > 2, the so-called strong collapse occurs, where the focusing 
happens with a power law dependence [12,13]. Hence, the critical case of the nonlinear focusing phenomena in equation 
(9) shares some of the collapse characteristics of the present model. 

We now turn to the massive case (still for k = 0). In the latter, we can no longer solve the problem analytically. 
However, we can directly observe that if we adopt a solution of logarithmic dependence of the form of Eq. (8) in 
this case as well, the massive terms diverge much more slowly (at worst as (log(f* — t)) 2 /(t* — t) 2 , as opposed to 
the (t* — t)~ 4 divergence of the dominant (massless) terms). This signifies that asymptotic self-similarity will ensue 
from this case and the logarithmic dependence will eventually set in and dominate the asymptotic behavior leading 
to collapse. An example of this is shown in Fig. 1, where a numerical simulation of Eq. (7) clearly indicates the 
logarithmic focusing of the scalar field for as t — > 1.38. This asymptotic type of self-similarity occurs quite often in 
cosmology, as can be seen for example in [14] (and references therein). One final comment worth making about this 
case is that it appears as if the massive case is asymptoting to the results for the massless one. 
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FIG. 1. Logarithmic blowup of the scalar field as a function of time in the massive case (for m = 1), starting with 0(0) = 1, 
0(0) = 0. The blowup occurs for t ~ 1.38. 

We now turn to the case with k ^ 0. In this case also, the reduction that leads to an equation only for the scalar 
field can be performed. However, due to the more complex nature of the equations, this reduction no longer leads to 
a second order ODE but rather to a third order one. In particular, in this case the reduction (after differentiating (4) 
and substituting the result, as well as (2) and (4), in (3)) results in 

600 - 80 2 - 1Oto 2 00 - 60 4 + 6m 2 2 + 3m 2 2 2 - 2to 4 2 = (10) 

It is interesting to note then that the massless case once again shares the exact same, finite time blowup solutions 
of Eq. (8). One can then once again use the same argument for the massive case to identify such solutions as the 
dominant asymptotic behavior for m =fi 0, since these terms blowup as (t* — t)~ 4 in Eq. (10), while the rest of the 
terms diverge with a rate of (at most) 0(t* — t)~ 2 . Numerical integration of Eq. (10) for various initial conditions 
and various masses confirms the theoretical prediction of finite time collapse. It appears however that the mass of the 
scalar affects the time at which collapse will occur (t*) (see e.g., Fig. 2). 




FIG. 2. Logarithmic blowup of the scalar field as a function of time in the massive case (for m = 0.1), starting with 0(0) = 1, 
0(0) = 0.5, 0(0) = (top panel) and 0(0) = 1, 0(0) = -0.5, 0(0) = (middle panel). The blowup occurs for t « 3 in these 
cases. A case with a larger mass (m = 0.6) is shown in the bottom panel. In this case, the blowup occurs at later times (in 
this case for t « 20) . 

It is also worth noting that in the context of Eq. (10), the steady state = appears to be linearly stable since 
for small perturbations ~ eexp(ivt), collecting the leading order behavior (0(e 2 )), one obtains the algebraic stability 
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equation of the form: 

(U? + TO 2 ) 2 = 0. (11) 

Eq. (11) has double imaginary roots: 

u) = ±im (12) 

which, in turn, denote marginal stability. Notice that the same situation occurs in the critical case of Eq. (9). However, 
the nonlinear term of the form — 60 4 asymptotically dominates and gives rise to a nonlinearity-induced instability 
leading to logarithmic blowup. This (blowup), in fact, is to be expected since the incompleteness of geodesies in 
3 + 1-dimensional spacctimcs [15], established using general techniques of differential topology in [16,17], has been 
shown under quite general conditions [18] to lead to divergence of physically observable properties. 

In conclusion, in this work, we have examined the coupling of 3 + 1-dimcnsional gravity to a scalar field satisfying a 
linear Klein-Gordon equation. We have found that it is possible to decouple the equation for the massive scalar from 
the one for the scale factor of the FRW metric used in this work at the "expense" of obtaining a nonlinear ordinary 
differential equation for the scalar. The "memory" of the coupling to gravity has been encapsulated in the nonlinear 
nature of the resulting equation. Closed form solutions of the resulting equation can be obtained in the massless 
case and exhibit logarithmic divergence of the scalar field as a function of time (and power-law vanishing of the scale 
factor). It is then observed that these solutions persist as dominant asymptotic behavior in the case where m =/= 0. 
These results are corroborated by numerical integration of the nonlinear ODEs in the massive case. Even though 
slightly different methods have been used for the cases where the metric is flat (k = 0) and when k ^ 0, the same 
principal conclusions have been drawn in all cases. 

This phenomenon of decoupling persists even in the case of more general, anisotropic Bianchi models [19]. In all 
these cases, special solutions to the ensuing decoupled non-linear ODE exist and they also exhibit blowup behavior. 

Finally, we return to the analogy of the logarithmic focusing of the solutions in this simple model (where spatial 
dependence was a priori fixed in the metric) with the logarithmic blowup in the critical case of a prototypical nonlinear 
partial differential equation, namely the nonlinear Schrodinger equation, that sustains focusing solutions. It would be 
particularly interesting to examine whether the inclusion of spatial dependence can lead to an equation of this type 
(through appropriate envelope wave expansions, given that NLS is the envelope equation for nonlinear wave equations 
of the Klein-Gordon type). In particular, if, in the presence of the spatial dependence, the reduction to a nonlinear 
partial differential equation for the scalar field provides a nonlinear KG equation, then it will be natural to expect 
that the reduction to NLS and the ensuing focusing solutions will carry through. If such a program succeeds, it may 
appear natural to consider the NLS as prototypical model for a "nonlinear quantum cosmology" . 
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